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I»    Introduction 

/'   If  a  liquid  (or  solid)  drop  is  siirrounded  by  a  gas  in 

which  the  pressure  of  the  drop's  vapor  is  below  the  equili- 
brium vapor  pressure,  the  vapor  diffuses  away  from  the  drop 
surface  and  the  drop  evaporates  (or  sublimates)  .   This 
evaporation,  which  maintains  the  vapor  pressure  near  the 
surface  at  the  equilibrium  value,  results  in  the  decay  of 
the  drop.   Conversely,  if  the  vapor  pressure  in  the  sur- 
rounding gas  exceeds  the  equilibrium  value,  vapor  diffuses 
toward  the  dr'op,  condenses  en  its  surface  and  the  drop 
grows.   The  same  description  applies  to  a  drop  dissolving 
or  grox-jing  in  a  liquid,  if  the  term  "vapor  pressure"  is 
replaced  by  "concentration  of  solute," 

We  wish  to  determine  the  drop  radius  and  the  vapor 
pressure  distribution  as  functions  of  time  when  a  single 
spherical  drop  evaporates  or  grows  in  an  infinite  motion - 
less  gas»   In  order  to  do  this,  we  shall  formulate  the 
physical  problem  mathematically  in  section  II.   This  prob- 
lem requires  that  the  diffusion  equation  for  the  vapor 
pressure  or  concentration  be  solved,  subject  to  c  onditions 
on  an  unknoi^m  moving  boundary,  vrhich  must  also  be  found. 
Previous  workers  have  approximated  the;  problem  either  by 
assuming  the  boundary  to  be  stationary  or  by  employing 
the  static  diffusion  equation.   By  treating  the  exact  prob- 
lem, we  find  that  the  first  of  these  approximations  is 
good  initially,  while  the  second  is  good  at  later  stages. 
Our  results  are  obtained  in  sections  III  and  compared  with 
these  approximations  in  section  IV, 

Although  the  equilibrium  vapor  pressure  or  concentration 
depends  upon  the  radius  and  the  temperature  of  thu  drop, 
we  assiwie  it  to  be  constant.   Thus  our  mathematical  problem 
aan  also  be  Interpreted  as  one  of  phase  change  produced  by 
heat  flov;,  provided  we  interpret  concentration  as  tempera- 
ture, diffusion  coefficient  as  heat  conductivity,  equili- 
brium concentration  as  transition  temperattire,  etc.   For 
example,  the  melting  of  ice  in  water,  or  of  any  solid  in 
a  liquid,  is  a  problem  of  this  type. 
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In  a  forthcoming  paper  i^ao   intend  to  present  the  extension 
of  the  present  methods  to  a  collocti  on  of  drops,  and  in  an- 
other, to  give  the  proof  of  the  existence  and  lonlqueness  of 
the  3D  lution  of  the  mathematical  problem  considered  herein. 

II o  Formulation 

VJe  consider  a  single  spherical  drop  of  radius  R(t)  at  time 
t,  which  has  radius  a  at  time  0,   If   c(r,t)  denotes  the  con- 
centration of  vapor  at  radial  distance  r  from  the  center  of 
the  drop  at  time  t,  then  we  assume  that  c  satisfies  the  dif- 
fusion equation 

(1)  °rr  +  f  °r  =  ^"^  °t   ' 

Here  D   is  the  coefficient  of  diffusion.   We  further  assume 
that,  at  the  surface  of  the  drop,  the  concentration  is  equal 
to  the  equilibrium  concentration  corresponding  to  the  tempera- 
ture and  curvatuxe  of  the  drop,  which  we  denote  by  g(t,R)  • 
We  assume  g(t,R)   is  a  known  function,  and  will  usually  as- 
sume it  to  be  constant.   Thus  we  have 

(2)  c(R(t),t)  =  g(t,R)   . 

The  rate  of  grov/th  or  decay  of  the  drop  depends  upon  the 
rate  at  which  vapor  condenses  on,  or  evaporates  from  it.   If 
p  denotes  the  density  of  the  liquid  in  the  drop,  then  conser- 
vation of  mass  yields 

Finally  we  assume  that  the  initial  concentration  is  uniform 
and  has  the  value   c„  .   Formally  we  have 

ik)  c(r,0)  =  Cq  . 

We  list  also  the  initial  condition  on  R(t),  which  was  stated 
above : 
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(5)      R(0)  =  a, 


The  problem  we  consider  is  that  of  determining  R(t)  and 
c(rpt)   for  subsequent  times,  and  in  particular,  of  finding  the 
time  required  for  a  droplet  to  disappear  completely  when  eva- 
poration occurs. 

To  this  end  we  first  introduce  new  dimensionless  variables 
as  follows: 

't  =  D  a'^t  c(F,t;  =  (gQ  -  CQ)"^[c(r,t)  -  c^] 

(6)  .^r  =  a'^r  h(t,R)  =  (g^-  c^)"^  [g(t,R)  -  c^] 


^R(t)  =  a-^R(t)  a   =  p"-^(gQ  -  Cq). 

In  these  equations  the  barred  quantities  are  the  new  dimension- 
less  variables,  g,  is  the  value  of   g(t,R)   at   t  =  0  ,   (1-5) 
now  become,  upon  introducing  (6),  the  following  (all  bars  have 
been  omitted)  : 

(7)  °rr  ^  §  °r  "  °t  '        ""  ="  ^^^^ 

(8)  c(R(t),t)  =  h(t,R)   , 

(9)  11-=^  c^(R(t),t)  , 

(10)  c(r,0)  =  0   ,  r  >  1 

(11)  R(0)  =  1   , 

In  the  case  of  condensation  (a  <  0)   R(t)   will  be  an  increas- 
ing function  of  t,   hence   R(t)  >  1  >  0  for  all  t,  and  the  region 
in  which  c(r,tj   is  desired  will  not  contain  the  line   r  =  0  • 
In  the  case  of  evaporation,  (a  >  0  ).  it  is  expected  that  R 
will  vanish  at  some  finite  time   t„  •   We  show  that  in  this  case 
an  additional  boundary  condition, 

(12)  0^(0, t)  =  0   ,       ^  >  ^0 

must  be  satisfied.   Indeed,  since  the  line   r  =  0  is  included  in 
the  region  of  definition  of  c(r,t),  this  function  can  be  extended 
for  r  negative.   Since  our  problem  has  spherical  symmetry,  c(r,t) 


must  be  an  even  function  of  r,  i.e^ 

c(r,t)  =  c(-  r,t)  . 

Continuity  of  c   and  of  c   in  the  region  above  the  evaporation 
curve  implies  then  the  condition  (12)  . 

Ill,   Solution 

Let  us  first  consider  the  one-dimensional  analogue  of  the 

2 
present  problem.   The  only  difference  is  that  tne  terra  —  c 

is  emitted  from  (1)  and  (7)".   If  h(t,R)=  1,  this  problem 
can  be  solved  exactly  and  has  the  solution 

,00    2   n"^  C^         2 


e"^  dx 


(13)     c(r,t)  =■/   e-^  dx 

i-'^o 
(111)     R(t)  =  1  +  2z^  \/t" 

where  the  quantity  z  is  defined  by 

(15)     z  =  (r  -  l)/2  it-   , 

and  z   is  the  unique  solution  of 

^2  ^00 
(16) 

If  a  <  0,   z   is  positive  and  the  droplet  grows.   If  a  >  0, 
then   z   <  0,  and  the  drop  evaporates,  in  which  case  solution 
(13-1[|.)  applies  until  the  drop  is  gone.   The  time  required  for 
complete  evaporation  is  obtained  from  (lU)  ^y  setting  R  =  0  . 
We  find 

(17)  t^  =  (2Zo)"^   • 

For  t  >  t  ,  condition  (12)  must  be  also  satisfied  and  this  is 
o' 

achieved  by  adding  to   c(rpt)  equation  (13),  ^  compensating 
part  u(r,t)  defined  as  the  solution  of 

(18)  u   =  u^^ 
^       '  rr    tt 

(19')      u(r,t^)  =  0,    r >  o 


V 
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(20)    u^CO.t)  =  -  0  ^   Jf    e-'d.  )-l  1    e-^/'^*,  t>t, 


Tf 


The  compensating  part  can  be  easily  constructed,  but  since  it 
is  felt  only  for   t  >  t  ,  it  will  not  affect  the  determination 
of  the  evaporation  curve. 

For  small   a,  the  approximate  solution  of  (l6)  is 

(21)       2,  =  -^  . 

Thus  the  time  required  for  evaporation,  if  a  is  small,  is  ap- 
proximately  (for  the  one -dimensional  case) 

2 
(22) 


The  last  expression  is  in  dimensional  unite. 

To  solve  fhe  spherical  problem,  we  employ  a  solution  of  (7) 
and  (10)  which  is  similar  to  (12),  but  involves  an  arbitrary 
function  f  which  is  to  be  determined.   This  solution  is 

itCO  p 

(23)      c(r,t)=i/   f  [(z-x)  VF]  e"^  dx   , 


z  being  defined  by  (l5)  .   (The  function  f  is  closely  related 
to  the  initial  value  of   c(r,t)   for  r  <  1;  precisely. 


f(^)  -^^   , 


r  <  1  .) 


In  the  case  of  evaporation,  the  additional  boundary  condition 
(12)  must  bo  satisfied  after  the  drop  disappears.   To  satis.'^y 
this  condition  we  need  merely  add  to  (23)  a  similar  expression 
with  -r   instead  of  r  after  the  drop  disappears.   That  this  is 
so  is  shovm  in  the  pq)  et  referred  to  above,  viiich  contains  the 
existence  and  uniqueness  proof.   This  compensating  part  plays 
no  role,  hovjovor,  in  the  determination  of  the  evaporation  curve. 

To  determine   R(t)  and  f('tr),  we  assume  that  R(t)  may  be  ex- 
panded in  a  formal  powej?  series  in  Vt  and  f(r)  in  a  formal 
power  series  in  tr  ,  i.e.  , 
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CD 


-n/2 


00 


n=l 


The  expression  "formal  power  series"  is  used  since  it  will  not 
be  established  whether  any  of  the  series  (2L|.)  and  (2^)      has  a 
positive  radius  of  convergence.   However,  it  will  be  demon- 
strated that  our  procediire  leads  to  results  satisfactory  for 
most  practical  purposes,  regardless  of  whether  the  above  series 
converge  or  do  not  converge.   In  any  case,  the  expression  (Slj.) 
for  R(t)  can  not  converge  up  to   t  =  t  ,  the  time  at  which  the 

drop  disappears,  and  therefore  could  not  be  used  for  determl- 

2 
nation   of  t   .   This  is  evidenced  by  the  fact  that  R  (t) 

behaves  like  a  linear  function  of  t  near   t  =  t   as  is  shown 

o 

in  the  paper  referred  to  in  the  Introduction,  (this  is  also 
confirmed  by  experimental  results)  and  a  series  of  th©  same 
type  as  in  (2[j.)  cannot  exhibit  such  a  behavior. 

By  means  of  (25),  equation  (23)  can  be  rewritten  as 


(26)   rc(r,t)  =  A 


.'/ 


CO 

e 


e"^  dx  + 


00 

rPl 


^n^ 


n/2 


(z-x)  e^  dx) 


a  and  A  are  now  determined  by  inserting  (26)  in 
n      n  ''  o    ^       ' 


The  numbers 

(8)  and  (9)  and  by  identifying  the  coefficients  of  various 

powers  of  yT  .   Thus  we  find,  when  h(t,R)  =  1 


(^7) 


/  ^1 


/ 


=  z 


a. 


A^  = 


Ui  = 


g  +  2(a+l)  ^o 

a+2(a+l)(l+z^'^) 

I'Oo     2    , 
(i    e"^  dx)"-^ 


I 


o 

a-a2 

a  +  1 


etc,    z      being    defined  by  equation    (16).    For   a   small,   we   have 
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(28) 


a. 


a 
2 


a 


)  ^3     "  3  yr  ' 


a 


V 


h 


(10+3Tt)a^ 


TStT 


""^^    \7¥=  ' 


A-,    =  3a    , 


a.-,   = 


-  2a 


lj.a' 


3 


As   already   stated,    the    series    {211)   ^^^  RCt)    can  at  best 
yield   information  about   the   local  behavior      of  R(t)   for   t   sraall. 
To   obtain  an  approximate   result   over   the  full   range    of  variation 
oft,      0<t<t,we   now   define    the   N'th  approximation 
Ojj(r,t)    to   c(r,t)      by 


(29)        rcj^(r,t)   =  A^( 


^00         2  N  /o 

e-^  dx  +  Ha  t^/2 

z  J 


CD 


(z-x)^e"^  dx)    , 


with  the  A  previously  computod( 


'jn 


The  N'th  approximation  RT,T(t)  to  R(t)  is  defined  by  requlimg 

J ''J 

that  the  mixed  condition 


(30)   [rc(r,t)]^  =  c(r,t)  +  rc^(r,t)  =  l  +  i  R  ||  , 


a 


r   =  R(t)  , 


derived  from  (8)  and  (9),  be  satisfied  exactly  by   c  =  c^^  and 
R  =  R. 


'N 


^N 


The  equation  determining  Rvr(t)  is  thus 


(31)  ^-k^t  (4^^))  =  ^o(-2Tr  "' "  i  t  -n^ 


n-1  r^      n-1  -x^ 

(z-x)^  -^e  ^  dx) 

z 


where 


z  = 


R^(t)  -  1 


2  VT 


and  Rjj(O)  =  1  . 


C„  and  R„  so  defined  satisfy  eqiations  (7),  (10),  and  (11),  and 
the  Slim  of  (8)  and  (9)  exactlyo   Furthermore  (8)  and  (9)  are  satis- 
fled  with  an  accuracy  increasing  with  N  for  small  t,  and  for  every 

2 
N,   R,T('t)  is  linear  in  t  near   t  =  t   , 

Equation  (31)  has  been  integiated  numerically  for  N  =  1,2,3 
and  a  =  2  X  10  "^  ,  corresponding  to  the  case  of  water  drop  in  air 
at  temperature  of  73°P  ,  assuming  the  vapor  pressure  at  infinity 
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to  be  ?xero.   The  results  listed  in  the  table  below  show  an  excel- 
lent rate  of  convergence. 

Since  z(0)  =  z  '^' ^  ,    and   z(to^  ~  "  \ *  ^^   ^^   expected 

that  z  is  a  small  n\iinber  for  aL  1  t  and  piay  be  replaced  in  (31)  by 
zero,  at  least  for  N  small.  The  r  esiO. ting  equation  for  R„  which 
approximates  R^(t)   is 

and  yields 

(33)   \  =  l-2at  +  2aA^(.  VI  +  J^^-lf   2H72  Hf  )V^  ^   * 

Comparison  with  results  of  nvimorical  integration  of  the  more  accu- 
rate equation  (31)  show?  again  a  good  agreement. 

From  (33  )>  ^e  easily  deduce  that  in  the  three-dimensional  case 

i3k)  t^  =ZS  &  -  3^  V^^  °(^))' 

and  that  to  a  good  approximation, 

(35)     R^(t)  =  1  -  2at  -  A  a  Vr^   for  0  <  t  <   to   . 

The  table  below  lists  the  values  of  R  (t)  computed  from  the 
series  (2L|.),  R^(t),  Rp(t),  and  Rr(  t  ),R  ( t)  tof  our  significant  figures, 

for  a  =  2  X  10,  .   To  this  approximation,  RT(t)  =  Rp('t),  and 
R^(t)  =  R2(t)  . 


•^ 


*s, 
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t        iR^icq, 2l\.) 


R- 


,5   X   10^ 

1  . 

2  ' 

l& 

8 
10 

12;   • 

Ik 
16- 

18. 
20. 

22. 

2[pcl0^ 
t_ 


1 
.9791 
.9590 

.9197 
.81^36 
.7712 
.7021 
.6362 

.5735 
.51U2 

.kSQi 
.i|o53 
.3557 
.3091+ 
,2662 


1 

.9790 

.9586 

.9180 

.8370 

.7563 

.6757 
f 

.5952 

.511+8 

.l;3l|4 
.35i]-0 
.2736 
.1932 
.1126 
.0325 
2,i|8lxlO 


R2 


«2 


R 


1 


-2 
«2 


h 


1 
.9790 
.9586 
.9181 
.8370 
.7562 
.6755 
.591+9 
.511+3 
.1+337 
.3531 
.2725 
.1919 
.1113 
.0307 


1 

.9790 
.9586 

I  .9180 

! 

I  .8371 

I 

i  .7565 

! 

I  .6760 

I  .5955 

I  .5151 

!  .1+3U7 

:  .351+3 

I  .2739 

1 

I 

i  .1936 

i 

i  .1133 

I 

i  .0330 


1 
.9790 

.9586 
.9181 

.8373 
.7568 
.676U 
.5961 
.5159 
.1+357 
.3555 
.2751+ 
.1953 
.1153 
.0352 


^1+ 


.h 


Formula    (3l+)   yields 


2.1+76x10^  ! 2. 1182x10"^ I  2.14-38x10^  | 

J+ 


for   t      the  value      2,1+88  x  10^   , 
o 
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rv  Comparison  With  Previous  Approximate  Solutions 

An  approximatlc^n,  expected  to  be  valid  initially,  can  be 
obtained  by  assuming  that  R(t)  s  1   and  solving  (7),  (S),  (10) 
for  c(r,t)  .   From  this  solution  c  (l,t)  can  be  computed  and 
used  in  (9)  to  determine  an  improved  R(t)  .  This  yields  the  be- 
ginning of  a  solution  in  powers  of  a  • 

We  easily  find  that 

p   / 
(36 


)   c(r,t)  .  ^,£   e-  h(t.  (^f  ,i;cl.  , 

(37)    E{t)    =   1    ^  ^Xjt      -  2^    I     h^  (t,l)dt  -  a 
If  h(t,l)  =  1,   this  becomes 


t  --t 

h(t,l)dt  . 
o 


(38)   R(t)  =  1  -  ^  yT  -  at 


• 


Comparison  i-jith  (2L(.)  and  (28)  shows  that  we  have  here  the  be- 
ginning of  the  asymptotic  expansion  for  R(t)   near  t  =  0  , 
provided  a  is  small. 

Another  approximation,  the  so-called  quasi-static  approxi- 
mation is  expected  to  be  valid  in  the  later  stages  of  the  process. 
It  is  obtained  by  first  neglecting  the  term  c,  in  (7)  and  the 
Initial  condition  (10).   Then,  from  (7),  (8)  and  from  the  as- 
sumi^tion  th^.t  c  tends  to  zero  as  r  becomes  infinite  •  (which  r- 
places  (10)),  we  find 

(39)    c(r,t)  =  h(t,R)R(t)   ^ 

Using  (9),  we  find 

(1^.0)   R^(t)  =  R^(t)  -  2a    h(t,R(t))dt,    -T  arbitrary  . 

This  is  an  integral  equation  for  R(t)  if  h  depends  upon  R,  but 
is  an  explicit  solution  if  h  is  independent  of  R,   VJith  R(t)  so 
determined,  c. (r,t)  may  be  computed  and  used  again  in  (7)  to 
obtain  an  improved  c(r,t)  •  V.'hen  h  =  1,   (i(-0)  yields 
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(i;l)   R^(t)  =  R^ra  -  2a(t  -  Z) 
=  1   -  2at   +   b   , 

where  b  =  R^  (tr)  +  2a 'T-  1  ,   If  T  =  0,  then  b  =  0  by  con- 
dition (11),  although  this  value  of  b  cannot  be  used  since 
the  quasi**static  solution  does  not  apply  initially,   i^ctual- 
ly,  comparison  with  (35)    shows  that  a  good  agreement  in  the 

later  stages  is  found  if  b  is  chosen  as   -2  (/-j^.  Physically 
this  means  that  the  drop  first  evaporates  quite  rapidly 
until  the  static  concentration  distribution  is  established 
and  then  behaves  quasi-statlcally,  starting  with  a  smaller 
radius. 

From  the  fact  that  the  first  approximation  above 
agrees  with  experimental  data  on  the  initial  growth  (decay) 
of  particles,  and  also  agrees  with  our  solution,  we  conclude 
that  our  solution  also  agrees  with  the  d^ta  initially.   Si- 
milarly, since  at  later  times  the  quasi-static  approxi- 
mation agrees  with  evaporation  data  and  with  our  solution, 
we  conclude  that  our  solution  agrees  with  these  data.   Thus 
we  expect  that  our  solution  will  agree  with  experiment  over 
the  entire  time  range  in  either  evaporation  or  growth. 
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